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ABSTRACT 

GALA is a freely distributed Fortran code to derive automatically the atmospheric parameters (temperature, 
gravity, microturbulent velocity and overall metallicity) and abundances for individual species of stellar spectra 
using the classical method based on the equivalent widths of metallic lines. The abundances of individual 
spectral lines are derived by using the WIDTH9 code developed by R. L. Kurucz. GALA is designed to obtain 
the best model atmosphere, by optimizing temperature, surface gravity, microturbulent velocity and metallicity, 
after rejecting the discrepant lines. Finally, it computes accurate internal errors for each atmospheric parameter 
and abundance. The code permits to obtain chemical abundances and atmospheric parameters for large stellar 
samples in a very short time, thus making GALA an useful tool in the epoch of the multi-object spectrographs 
and large surveys. An extensive set of tests with both synthetic and observed spectra is performed and discussed 
to explore the capabilities and robustness of the code. 
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1. INTRODUCTION 

The last decade has seen a significant improvement in 
the study of the chemical composition of the stellar pop- 
ulations (in our Galaxy and its satellites), thanks to the 
8-10 meters class telescopes, coupled with the design of 
several multi-object mid/high-resolution spectrographs, e.g. 
FLAMES mounted at the Very Large Telescope, AAOmega 
at the Anglo- Austi-alian Telescope, DEIMOS at the Keck Ob- 
servatory and HYDRA at the Blanco Telescope of the Cerro 
Tololo Inter- American Observatory. These instruments have 
allowed to enlarge the statistical significance of the acquired 
stellar spectra, but also they have required a relevant effort to 
manage such large databases. 

The next advent of new surveys aimed to collect huge sam- 
ples of mid-high resolution stellar spectra, as for instance 
the European Space Agency Gaia mission, the Gaia-ESO 
Survey committed b y the European Southern Observatory 
(lGilmoreetal.1 120121). the APOGEE Survey at the Apache 
Point Observatory (lAllende Prieto et al.l2008l). and the RAVE 
Surve y at the Anglo- Australian Observatory (ISteinmetz et alj 
120061) . will make available in real time to the astronomi- 
cal community an enormous volume of data. Other spec- 
troscopic surv eys have been alread y performed, like B RAVA 
( iKunderetaij|2012l) and ARGOS ( Freeman et aUlMSi) . both 
dedicated to the study of the Galactic Bulge. Also, other 
multi-object spectrographs are planned or p roposed for the 
next years, i.e. HERMES (Barden et al."2010) at the Anglo- 
AustraHan Observatory, 4MOST (de Jong 2011) at the New 
Technology Telescope and MOONS (Ciras uolo et al.ll201 II) 
at the Very Large Telescope. This perspective, coupled with 
the huge amount of high quality spectra available in the main 
on-line archives (and not yet totally analysed) clarifies the ur- 
gency to develop automatic tools able to rapidly and reliably 
manage such samples of spectra. 

In the last years several codes aimed at the automatic 

' Based on observations collected at the ESO-VLT under programs 
65.L-0165, 165.L-0263, 073.D-0211, 080.D-0368, 083.D-0208 and 
266.D-5655 and on data available in the ELODIE archive. This research 
has also made use of the SIMBAD database, operated at CDS, Strasbourg, 
France. 



measurements of the chemical abundances have been al- 
ready developed. They are mainly based on the com- 
parison between the observed spectrum and grids of syn- 
thetic spect ra, for instance ABB O (Bonifacio & Caffau 2003), 
MATISSE (Recio-Blanco, Bijaoui & de Laverny 2006), SME 
(Valenti&Piskunov 1996), SPADES (Posbicetal. 2012), 
MyGIsFOS (Sb ordone et al.i i2010). In particular, in these 
codes the main effort has been devoted to robustly determine 
the atmospheric parameters (and hence the elemental abun- 
dances) for low (<50) signal-to-noise (SNR) spectra and gen- 
erally to develop an algorithm able to accurately treat different 
kind of stars (in terms of metallicity and stellar parameters). 

In this paper we present and discuss a new code (named 
GALA) specifically designed for automatically determining 
the atmospheric parameters by using the observed equiv- 
alent widths (EWs) of metallic lines in stellar spectra, 
at variance with the majority of the available automated 
codes. GALA is a tool developed within Cosmic-Lab, a 
5-years project funded by the European Research Coun- 
cil and it is freely available at the website of the project 
^http://www.cosmic-lab.eu/Cosmic-Lab/Products.html 

The paper is structured as follows: Section 2 discusses the 
outline of the classical method to derive the main parameters 
and the interplay occurring among them; Section 3 describes 
the algorithm; Section 4 describes the identification and the 
rejection of the outliers and Section 5 discusses other aspects 
of the code. Section 6 provides a complete description of the 
uncertainties in the calculations. Finally, Sections 7, 8 and 9 
discuss a number of tests performed to check the stability and 
the performances of GALA. 

2. THE METHOD 

The main advantage of inferring the stellar atmospheric pa- 
rameters from the EWs is its reproducibility: any researcher 
can directly compare its own results about a given star with 
other analysis based on the same approach. This allows to 
distinguish between discrepancies due to the method (i.e. the 
measured EWs) and those due to the physical assumptions of 
the analysis (model atmospheres, atomic data...). On the other 
side, one of the most critical aspects of this method is the 
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particular accuracy needed in the definition of the lineHst, by 
excluding blended lines. In fact, the codes developed to calcu- 
late the abundance from the measured EW compare the latter 
with the theoretical strength of the line, changing the abun- 
dance until the observed and theoretical EW match within 
a convergence range. The theoretical line profile is usually 
calculated including the continuum opacity sources but ne- 

f' ;cting the contribution of the neighboring lines (see lCasteliU 
05a, for details), hence the spectral lines to be analysed 
with this technique are to be checked accurately against blend- 
ing (a practice not always performed). Otherwise, the use of 
synthetic spectra allows to also use blended features (and in 
principle to exploit the information derived from all the pix- 
els), but it is more expensive in terms of computing time, be- 
cause large grids of spectra must to be computed, at differ- 
ent parameters and with different chemical compositions, and 
each change in the atomic data leads to a recomputation of the 
synthetic spectra. 

2.1. The classical spectroscopic method 

The main parameters that define the model atmosphere, 
namely the effective temperature {T^a), the surface gravity 
(log g), the microturbulent velocity (vt) and the overall metal- 
Ucity ([M/Hfl) are constrained through four constraints: 

1. Temperature: the best value of Teff is derived by im- 
posing the so-called excitation equilibrium, requiring 
that there is no correlation between the abundance and 
the excitation potential x of the neutral iron lines. The 
number of electrons populating each energy level is 
basically a function of Teff, according to the Boltzmann 
equation. If we assume a wrong Teff in the analysis of 
a given stellar spectrum, we need different abundances 
for matching the observed profile of transitions with 
different x- For instance, the use of a value of Teft too 
large will lead to under-populate the lower energy lev- 
els, thus the predicted line profile for low-x transitions 
will be too shallow and a higher abundance will be 
needed to match the line profile. On the other hand, a 
wrong (too low) T^a will lead to a deeper line profile 
for the low x transitions. 

For this reason, a wrong, too large value of T^a will 
introduce an antic orrelation between abundances and 
X, and in the same way, a positive correlation is 
expected in the case of the adoption of a value of Teff 
too small; 

2. Surface gravity: the best value of log g is derived with 
the so-called ionization equilibrium method, requiring 
that for a given species, the same abundance (within 
the uncertainties) has been obtained from lines of 
two ionization states (typically, neutral and singly 
ionized lines). Being the gravity a direct measure of 
the pressure of the photosphere, variations of log g 
lead to variations of the ionized lines (very sensitive 
to the electronic pressure), while the neutral lines are 
basically insensitive to this parameter. This method 
assumes implicitly that the energy levels of a given 
species are populated according to the Boltzmann and 
Saha equations (thus, under Local Thermodynamical 

- We adopted the classical bracket notation where [X/H]= A(X)i,or-A(X)Q, 
where A(X)= logNA/Nn + 12. 



Equilibrium (LTE) conditions). Possible departures 
from this assumption (especially critical for metal-poor 
and/or low-gravity stars) could alter the derived gravity 
when it is derived from the ionization balance, because 
non-LTE effects affect mainly the neutral lines (even if 
the precise magnitude of the departures from LTE for 
the iron lines is still matter of debate). As s anity c heck, 
following the suggestion by lEdvardssonl (119881) . sur- 
face gravities determined from the ionization equilibria 
have to be "checked - when possible - with gravities 
determined from the wings of pressure broadened 
metal lines"; 



3. Microturbulent velocity: Vt is computed by requiring 
that there is no correlation b etween the ir on abundance 
and the line strength (see iMucciareilH |201 U for a 
discussion about different approaches). The microtur- 
bulent velocity affects mainly the moderate/strong lines 
located along the flat regime of the curve of growth, 
while the lines along the linear part of the curve of 
growth are mainly sensitive to the abundance instead 
of the velocity fields. The necessity to introduce the 
microturbulent velocity as an additional broadening 
(added in quadrature to the Doppler broadening) arises 
from the fact that the non-thermal motions (basically 
due to the onset of the convection in the photosphere) 
are generally not well described by the 1 -dimensional, 
static model atmospheres. Citing Kurucz (2005), 
microturbulent velocity is a parameter that is generally 
not considered physically except in the Sun, because in 
the Sun the velocity fields can be derived as a function 
of the optical depth through t he analysis of the intensity 
spectrum (as performed by Fontenla. Avrett & Loeserl 
|l993). For the other stars v't represents only a corrective 
factor that minimizes the line-to-line scatter for a given 
species and it compensates (at least partially) to the in- 
complete description of the convection as implemented 
in the 1 -dimensional model atmospheres; 



4. Metallicity: [M/H] is chosen according to the average 
iron content of the star, assuming [Fe/H] as a proxy of 
the overall metallicity. Generally, [Fe/H] is adopted as a 
good proxy of the metallicity because of its large num- 
ber of available lines, but it does not indicate necessar- 
ily the overall metallicity of the studied star. In fact, 
iron is generally not the most abundant element in the 
stars, while elements as C, N and O would be the best 
tracers of the stellar metallicity (but they are difficult to 
measure). 

Because of its statistical nature, the spectroscopic optimiza- 
tion of all the parameters simultaneously can be performed 
only if we have a sufficient number of Fe lines, distributed 
in a large range of EW and x ^nd in two levels of ioniza- 
tion. Alternatively, Teff and log g can be inferred from the 
photometry (for instance with the isochrone-fitting technique 
or employing empirical or theoretical Tgff-color relations) or 
by fitting the wings of damped lines (as the hydrogen Balmer 
lines or the Mg b triplet) sensitive to the parameters, and only 
Vt needs to be tuned spectroscopically (following the approach 
described above). Note that some authors consider the method 
to derive the parameters based on these constraints only as 
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sanity checks performed a posteriori on the photometric pa- 
rameters, while other authors rely on these constraints to infer 
the best parameters. 

2.2. The interplay among the parameters 

In light of the method described above, it is worth to bear 
in mind that the atmospheric parameters are correlated with 
each other In fact, the strongest lines are typically those with 
low x- Fig- [H plots all the transitions in the range A= 4000- 
8000 A and with x <10 eV in the Kurucz/Castelli databasfl 
in the plane x vs log(§/)-6'xQ. As mentioned in Section lTTl 
Teff and x strictly linked, and there is also a connection 
between Vt and the line strength. Hence, the statistical cor- 
relation between x ^nd the line strength leads to a correlation 
between T^a and Vt. Thus, a variation of Teff implies a variation 
of Vt. Also, variations of T^a and Vt will change differently the 
abundances derived from different levels of ionization (hence, 
the gravity). 
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Fig. 1 . — Behaviour of the hne strength (computed assuming reff=4500 K) 
as a function of the excitation potential x for the transitions available in 
the Kurucz/Castelli linelist in the range of wavelength A=4000-8000 A and 
with X <10 eV. 

Let us consider an ATLAS9 model atmosphere computed 
with 7; tt=4500 K, log g=1.5, v,=2 km/s and [M/H]=-1.0 dex 
(ICaste lli & Kumcz 2004) and a set of neutral and singly ion- 
ized iron lines (predicted to be unblended through the inspec- 
tion of a synthetic spectrum calculated with the same param- 
eters and convoluted at a spectral resolution of 45000). The 
EWs of these transitions are computed by integrating the the- 
oretical line profile throu gh the WID subroutine implemented 
in the WIDTH9 code (C astellil [2005a') . This means that each 
of these EWs will provide exactly [El/H]=-1.0 dex when they 
are analysed by using the model atmosphere described above. 

' http://wwwuseroat.ts.astro.it/castelU/linelists.html 
* The term log(g/) — 9x is used as theoretical proxy of the line strength, 
where log(gf) is the oscillator strength and 9= 5040/reff eV^^ . 



The analysis of these lines (adopting always the same set 
of EWs) is repeated investigating a regular grid of the atmo- 
spheric parameters, namely T^ff= 3600-5400 K, log g= 0.5- 
2.5, v,= 1.0-3.0 km/s, steps of 5%fs= 200 K, 5logg= 0.2, 
5v,= 0.5 km/s and assuming for all the models [M/H]= -1.0 
dex. Fig. 12] |3] and |4] display the quite complex interplay oc- 
curring among the atmospheric parameters. 

• Fig. |2] shows the behaviour of the slope of the 
A(Fe)-x relation for the above sample of lines as a 
function of Teft, keeping gravity fixed, but varying the 
microturbulent velocity. The thick grey line connects 
points calculated with the original Vt of the model. The 
global trend is basically linear, at least if we consider a 
range of ±1000 K around the original Tgff. The inset 
panel shows the behaviour of Tgff for which 5^ is zero 
(thus, the best Tgff) as a function of Vt and considering 
different gravities. The derived best temperature 
increases with increasing Vt (at fixed log g); the gravity 
has only a second-order effect and it does not change 
the general behaviour of the best Teff as function of Vt. 
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Fig. 2. — Main panel: behaviour of the slope 5^;. of the A(Fe)-x relation as 
a function of Teff assuming log g= 1 .5 and for different values of w (dashed 
curves). The thick grey curve represents the behaviour computed for the orig- 
inal microturbulent velocity (v,= 2.0 km/s). Horizontal dashed line is the zero 
value (according to the excitation equilibrium). The inset panel shows the be- 
haviour of the best value of T^ff (for which Sy^= 0) as a function of the used 
I't (empty squares) and for different gravities; black points are referred to the 
oiiginal gravity (log g=1.5). 

• Fig- E] summaryzes the behaviour of the slope Sewr 
of the A(Fe)-EWR relation as a function of Vt 
(where EWR indicates the reduced EW, defined as 
EWR=lg(£'W / A)), keeping gravity fixed at the original 
value but varying Tgff. For a given temperature, the 
slope decreases increasing Vt, with a behaviour that 
becomes less steep at Vt larger than the original value. 
The effect of the Teff is appreciable for Teft larger than 
the original value, while for lower Teff all the curve 
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are very similar to each other The inset shows the 
change of the best value of Vt (for which the slope 
of the A(Fe)-EWR relation is zero) as a function 
of Teff and for different log g. The observed trend 
is quite complex: basically, we note that the best 
value of Vt is very sensitive to T^ff, when the latter is 
overestimated with respect to the true temperature, but 
with a negligible dependence from gravity, while the 
behaviour is the opposite when Teff is under-estimated, 
with a degeneracy between T^ff and the best v, but a 
consistent dependence from gravity. 
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Fig. 3. — Main panel: behaviour of the slope Sewr of the A(Fe)-EWR re- 
lation a.s a function of v, assuming log g= 1.5 and for different T^d (dashed 
curves). The thick grey curve represents the behaviour computed for the orig- 
inal temperature (Tcff= 4500 K). Horizontal dashed line is the zero value. The 
inset panel shows the behaviour of the best value of the microturbulent ve- 
locity V, (for which Sewr= 0) as a function of the used T^ff (empty squares) 
and for different values of gravities; black points are referred to the original 
gravity (log g=1.5). 



• Fig- El shows the behaviour of the difference between 
A(Fe 1) and A(Fe 11) as a function of log g, assuming 
Vt= 2 km/s and for different Teff . The general behaviour 
is linear and the iron difference increases considerably 
increasing temperature. The best value of gravity (see 
the inset in Fig.|4|i is highly sensitive to changes in Teff, 
with Slogg/STsff ~1 dex/300 K in the investigated case, 
but for a fixed T^ff turns out to be marginally sensitive to 
Vt (this is due to the fact that the Fe 11 lines are basically 
distributed in strength in a similar way to the Fe 1 lines). 

It is important to bear in mind that the these considerations 
are appropriate for the investigated case of a late-type star but 
the dependencies among the parameters can be different for 
different regimes of atmospheric parameters and/or metallic- 
ity. However, the example presented above demonstrates that 
an analytic approach to derive the best model atmosphere is 
discouraged because it needs to know the precise topography 
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Fig. 4. — Main panel: behaviour of the difference between A(Fe I) and 
A(Fe II) as a function of log g assuming Vt= 2.0 km/s and for different T^f[ 
(dashed curves). The thick grey curve represents the behaviour computed for 
the original temperature (reff= 4500 K). Horizontal dashed line is the zero 
value (according to the ionization equilibrium). The inset panel shows the 
behaviour of the best log g (for which A(Fe I)= (Fe II)) as a function of the 
used Teff (empty squares) and for different gravities; black points are referred 
to the original microturbulent velocity (vt= 2.0 km/s). 

of the parameters space and requires the inspection of a large 
number of model atmospheres. 

3. GALA 

GALA is a program written in standard Fortran 77, that uses 
the W1DTH9 co de developed by R. L. Kumcz in its Linux ver- 
sion jSbordone e t al. 2004) to derive the chemical abundances 
of single, unblended absorption lines starting from their mea- 
sured EWs. We used our own version of W1DTH9, modified 
in order to have a more flexible format for the input/output 
files with respect to the standard version of the code available 
in the website of F. Castelli, while the input physics and the 
method to derive the abundances are unchanged. 

GALA is specificaUy designed to 

1. choose the best model atmosphere by using the ob- 
served EWs of metallic lines; 

2. manage the input/output files of the WIDTH9 code; 

3. provide statistical and graphical tools to evaluate the 
quality of the final solution and the uncertainty of the 
derived parameters. 

GALA is designed to handle bofli ATLAS9 
(Caste lU & Kuruczl |2004) and MARCS (Gustafsso n et alJ 
2008) model atmospheres, that are the most popular em- 
ployed dataset of models. The current version has been 
compiled with the Intel Fortran Compiler (versions 11, 12 
and 13) and tested on the Leopard, Snow Leopard and Lion 
Mac OSX systems, and on the Ubuntu, Fedora and Mandriva 
Linux platforms. 
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3.1. Optimization parameters 

GALA has been designed to perform the classical chemical 
analysis based on the EWs in an automated way. The user 
can choose to perform a full spectroscopic optimization of the 
parameters or to optimize only some parameters, keeping the 
other parameters fixed to the specified input values. 

The algorithm optmizes one parameter at a time, checking 
continuously if the new value of a given parameter changes 
the validity of the previously optmized ones. 

For each atmospheric parameter X (corresponding to Teff, 
log g, Vt and [M/H]) we adopt a specific optimization parame- 
ter C(X), defined in a way that it turns out to be zero when the 
best value of the X parameter has been found. Hence, GALA 
varies X until a positive/negative pair of the C(X) is found, 
thus bracketing the zero value corresponding to the best value. 
Thus, the condition C(X)= identifies X=X as best value of 
the given parameter. Finally, the best solution converges to 
a set of parameters which verifies simultaneously the con- 
straints described in Section 12.11 When the atmospheric pa- 
rameters have been found, the abundances of all the elements 
for which EWs have been provided are derived. 

According to the literature, the adopted C(X) hav e bee n de- 
fined to parametrize the conditions listed in Section l2n 

1. the angular coefficient of the A(Fe)-x relation {S^} 
to constrain Teff; the lower panel of Fig. |5] shows the 
change of this slope for a set of theoretical EWs com- 
puted with the correct T^ff (grey points) and with tem- 
peratures varied by ±500 K (empty points). The vari- 
ation of Teff produces a change of the slope (but also a 
change in the y-intercept); 

2. the angular coefficient of the A(Fe)-EWR relation 
(Sewr) to constrain Vt. The upper panel of Fig.|5]shows 
the same set of theoretical EWs analysed with different 
Vt; 

3. the difference of the mean abundances obtained from 
Fe I and Fe II lines to constrain the gravity; 

4. the average Fe abundance to constrain the metallicity of 
the model. 

If the errors in the EW measurement ((Tew) are provided 
as an input, the slopes are computed taking into account the 
abundance uncertainties of the individual lines; the uncer- 
tainty on the iron abundance of a given line is estimated from 
the difference of the iron abundance computed for the input 
EW and for EW+aEw Q In the A(Fe)-x plane the uncertain- 
ties of X are reasonably assumed negligible, because the un- 
certainties of X are typically less than 0.01 eV, while in the 
A(Fe)-EWR plane the least square fit takes into account the 
uncertainties in b oth the axis (following the prescriptions by 
iPress et al.|[T992l) . 

The flexibility of GALA permits to simultaneously deal 
with stars of different spectral types. This is made by consid- 
ering that a large number of Fe I lines is generally available 
forF-G-K spectral types stars, whereas they are less numerous 
(or lacking) in O-B-A stars, for which a large number of Fe II 
lines is typically available. Also, in some spectral regions 
there is a large number of lines for other iron-peak elements 

' The uncertainties in A(Fe) ai'e assumed symmetiic with respect to the 
variations of the EW (ito-£iv); we checked that this assumption is correct at 
a level of ~0.01 dex. 
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Fig. 5. — Upper panel: behaviour of A(Fe I) as a function of the reduced 
equivalent widths for a set of theoretical EWs obtained from a model atmo- 
sphere computed with reff= 4500 K, log g= 2.0 and Vt= 2 km/s. The deiived 
abundances are obtained by adopting the correct value of I't (grey points) and 
two wrong values of V| (open circles). Lower panel: behaviour of A(Fe 1) as a 
function of the excitation potential for the same set of theoretical EWs. Grey 
points are the results obtained by analyzing the lines with the correct value 
of Tiff, while the empty circles are obtained by over/under-estimate T^ff by 
±500 K. 

(mainly Ni, Cr and Ti). For this reason, GALA is designed to 
optimize the parameters also using other lines instead of Fe I 
lines, by appropriately configuring the code. In the following, 
we will refer to the optimization made by using Fe I for T^ff 
and Vt, but our considerations are valid also for other elements 
with a sufficient number of lines. 

3.2. The main structure 

GALA is structured in three main Working-Blocks: 

(1) the Guess Working-Block is aimed at finding the guess 
atmospheric parameters in a fast way (this is especially useful 
in cases of large uncertainties or in lacking of first-guess 
value for the parameters); 

(2) the Analysis Working-Block which finds the best model 
atmosphere through a local minimization and starting from 
the guess parameters provided by the user or obtained through 
the previous block; 

(3) the Refinement Working-Block that refines the solution, 
starting from the atmospheric parameters obtained in the 
previous block. 

GALA can be flexibly configured to use different combina- 
tions of th e thr ee main Working-Blocks. We defer the reader 
to Section [84l for a discussion of the effects of the Working- 
Blocks. In the following we describe the algorithm of each 
block and the cases in which they are recommended. 

3.2.1. Guess Working-Block 

If the atmospheric parameters are poorly known, this 
Working-Block verifies them quickly by exploring the param- 
eters space in a coarse grid. Thus, it saves a large amount of 
time if the initial parameters are far away from the correct 
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solution. 

1 . As a first step, the abundances for each line are derived 
with the input parameters and some lines are labelled 
as outliers and excluded from the analysis (the criteria 
of the rejection are described in Sect. 0. The surviving 
lines will be used in this Working-Block and no other 
line rejection will be performed until convergence. 

2. The metallicity of the model is eventually readjusted 
according to the average iron abundance. 

3. 5^ is computed with the input T^ff and with a T^ff varied 
by +500K (if is positive) or -500 K (if is neg- 
ative). This procedure is repeated until a pair of posi- 
tive/negative values of is found, thus to bracket the 
Teff value for which 5'^= 0. The behaviour of 5^ as a 
function of T^ff is described with a linear relation, find- 
ing the value of T^ff for which 8^= 0. The description 
of this relation with a linear fit is legitimate as long as 
the employed Tgff range is relatively small (in this case 
500 K), because for larger range the behaviour of S^^ as 
a function of Teff could become non linear (mainly due 
to the interplay with the other parameters). 

4. The new value of Tgff is adopted to find a new value 
of Vt, following the same approach used for Teff and 
searching for a positive/negative pair of Sewr over a 
range of 0.5 km/s. 

5. Finally, a new value of log g is found, starting from 
the Teff and Vt derived above, by searching for a posi- 
tive/negative pair of A(Fe), over a range of 0.5 dex in 
gravity. 

The entire procedure (from the optimization of [M/H] to 
that of log g) is repeated for a number of iterations chosen by 
the user and finally a new set of input parameters is found. 
Generally three/four iterations are sufficient to find a good so- 
lution. The final solution is accurate enough to identify the 
neighborhood of the real solution in the parameters space but 
it could be unreliable since it needs to be checked for the 
covariances among the parameters (which is the task of the 
Analysis Working-Block. 

3.2.2. Analysis Working-Block 

This Working-Block performs a complete optimization 
starting from the input parameters provided by the users or 
from those obtained with the Guess Working-Block. This 
block is developed to find a robust solution under the assump- 
tion that the input values are reasonably close to the real so- 
lution. When good priors are available (for instance, in the 
case of stellar cluster stars) this block is sufficient to find the 
solution, without the use of the guess block. Otherwise, when 
the guess model is uncertain (for instance, in the case of field 
stars for whose reddening, distance and evaluative mass could 
be highly uncertain, or, generally, in the case of inaccurate 
photometry), the analysis block is recommended to be used 
after the guess block. Fig. |6] shows as a flow chart the main 
steps of this Working-Block. 

The following iterative procedure is performed: 

1 . The procedure starts by computing the abundances us- 
ing the guess parameters and performing a new line re- 
jection (independent from that of the previous block). 
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Fig. 6. — Flow diagram for the Analysis Working-Block of GALA (see 
Section 3.2 for details). 

At variance with the previous Working-Block, now the 
parameters are varied by little steps, configured by the 



2. The model metallicity is refined to match the average 
iron abundance. 

3. A new model with different Teff is computed, accord- 
ing to the sign of of the previous model (i.e., a 
negative slope indicates a overestimated Teff and vice 
versa). New models, varying only Teff, are computed 
until a pair of negative/positive is found. Thus, these 
two values of Teff identify the range of Teff where the 
slope is zero. Teff corresponding to the minimum \Sx \ is 
adopted. 

4. The same procedure is performed for Vt. If the final 
value of Vt is different from that used in the previous 
loop, GALA goes back to item (2), checking if the new 
value of Vt needs a change in [M/H] and Teff. Otherwise, 
the procedure moves on to the next loop. 

5. The surface gravity is varied until a positive/negative 
pairs of A{Fe) is found. If the output log g differs from 
the input value, GALA returns to (2) with the last ob- 
tained model atmosphere and the entire procedure is re- 
peated. 

6. When a model that satisfies simultaneously the four 
constraints is found, the procedure ends and the next 
star is analysed. 

We stress that the method employed in this Working-Block 
is very robust but it has the disadvantage of being slow if the 
guess parameters are far from the local solution. 

3.2.3. Refinement Working-Block 

This block allows to repeat the previous Working-Block 
using the solution obtained in the previous block as a starting 
point. A new rejection of the outliers is performed and the 
same approach of the Analysis Working-Block is used. This 
block can be useful to refine locally the solution when the 
first block is switched off. 

The main advantage of the Refinement Working-Block is 
that the new line-rejection is performed by using accurate at- 
mospheric parameters (because obtained from the Analysis 
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Working-Block). As will be discuss in Section |4] the line- 
rejection performed on abundance distributions obtained with 
wrong parameters can be risky, losing some useful lines. 

4. WEEDING OUT THE OUTLIERS 

The detection and the rejection of lines with discrepant 
abundances are crucial aspects of the procedure and require 
some additional discussion. Before ruling out a line from the 
line list, we need to understand the origin of the detected dis- 
crepancy. Basically, the main reasons for a discrepant abun- 
dance are: 

(1) inaccurate atomic data (i.e. oscillator strengths) that can 
under- or over-estimate the abundance; 

(2) unrecognized blends with other lines (providing systemat- 
ically overestimated abundances); 

(3) inaccurate EW measurement. 

The first two cases can be partially avoided with an effort dur- 
ing the definition of the adopted linelist, including only transi- 
tions with accurate log gf and checking each transition against 
blending, according to the atmospheric parameters and the 
spectral resolution. 
GALA rejects the lines according to the following criteria: 

1. lines weaker or stronger than the input EWR thresh- 
olds are rejected, in order to exclude either weak and/or 
strong lines. In fact, weak lines can be heavily affected 
by the noise, whereas strong lines can be too sensitive 
to Vt and/or they can have damping wings for which the 
fit with a Gaussian profile could be inappropriate, pro- 
viding a systematic under-estimate of the EW; 

2. lines whose uncertainty on the EW measurement is 
larger than an input threshold chosen by the user and 
expressed as a percentage. Note that not all the codes 
developed to measure EWs provide an estimate of the 
EW error, despite the importance of this quantity. For 
instance, among the publicl y available codes aimed to 
measure EWs, DAOSPEC (Stetson & Pancinol l2008h 
and EWDET (Ramirez et al. 2001) provide accurate 
uncertainty evaluations for each line, while SP ECTRE 
CFitzpatrick & Sneden 1987 ) and ARES ( Sousaet al.1 
120071) do^not include EW error calculations. For this 
reason GALA works even if aEw are not provided, al- 
though this affects the final solution accuracy, because 
all the transitions will be weighed equally, despite their 
different measurement quality; 

3. lines are rejected according to their distance from the 
best-fit lines computed in the A(Fe)-x and A(Fe)-EWR 
planes through a cr-rejection algorithm. A cr-rejection 
from the best-fit lines in the planes used for the op- 
timization is more robust with respect to a simple cr- 
rejection based on the abundance distribution. In the 
latter case, there is the risk to lose some lines important 
for the analysis, thus biassing the results. Fig|7]explains 
this aspect: we consider a synthetic spectrum of a gi- 
ant star (reff= 4500 K) and we measure the EWs after 
the injection of Poissonian noise in the spectrum in or- 
der to reproduce a reasonable good SNR (^30). Fig.|2] 
shows the distribution of the Fe I lines in the A(Fe)-x 
plane when the chemical analysis is performed by us- 
ing a wrong model atmosphere with reff= 5200 K (thus 
leading to an anticorrelation between A(Fe I) and x)- In 
the upper panel the outliers were rejected according to 
the median value of the abundance distribution, shown 



as gray solid line, while the two dashed lines mark ±.3a 
level and black points are the surviving lines. In the 
lower panel the rejection is performed according to the 
distance from the best-fit line (shown as solid line while 
the two dashed lines mark ±3cr level). It is evident that 
in the first case the majority of the discarded lines are 
those with low x (thus, the most sensitive to the T^ff 
changes), with the risk of introducing a bias in the Teff 
determination. On the other hand, the method of rejec- 
tion shown in the lower panel of Fig. |7] preserves the 
low-x lines, guaranteeing the correctness of the final 
solution. 
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Fig. 7. — Behaviour of the Fe I abundances as a function of the excitation 
potential for a synthetic spectrum computed assuming ri;ff= 4500 K but anal- 
ysed with a model atmosphere with T^ff= 5200 K. Black circles are the lines 
survived after the line-rejection procedure and the empty points the rejected 
lines. Upper panel shows the results by adopting a rejection based on the 
abundance distribution; the sohd Hne indicates the median abundance and the 
dashed lines mark ±3(7 level. The lower panel shows the results by using the 
procedure employed by GALA: solid line is the best-fit line and the dashed 
lines mark ±3(7 from the best-fit hne. 



An important point is that the outlier rejection in GALA 
is not performed independently in each iteration of the code, 
but only at the beginning of each Working-Block. This is es- 
pecially important, because it allows to use always the same 
sample of lines during the optimization process, avoiding the 
risk to introduce spurious trends in the behaviour of the given 
optimization parameter as a function of the corresponding at- 
mospheric parameters. In fact, the values of C(X) derived 
from two different sets of lines of the same spectrum, but for 
which an independent rejection of the outliers has been per- 
formed, cannot be directly compared to each other to derive 
X. In particular, this effect is magnified in cases of small num- 
ber of lines, where the impact of the lines rejection can be 
critical. 

5. MORE DETAILS 
5.L A comment about the gravity 
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The most difficult parameter to be constrained with the clas- 
sical spectroscopic method is the gravity. This because of the 
relatively small number of available Fe II lines, which can 
vary in the visual range from a handful of transitions up to 
^20, depending on the spectral region and/or the metallic- 
ity (for instance, some high-resolution spectra with a small 
wavelength coverage, as the GIRAFFE@VLT or the Hy- 
dra@BlancoTelescope spectra, can be totally lacking in Fe II 
lines). 

GALA is equipped with different options to optimize log g: 

(1) the normal optimization by using the difference between 
the average abundances from neutral and singly ionized iron 
lines (as described above); 

(2) the gravity is computed from the Stefan-Boltzmann equa- 
tion, by providing in input the term e = log(4GM7r(T/L), where 
G is the gravitational constant, a the Boltzmann constant and 
M and L are the mass and the luminosity of the star Thus, 
during the optimization process, the gravity is re-computed 
(as log g= e-41ogreff) in each iteration according to the new 
value of Teff; 

(3) the gravity is computed by assuming a quadratic relation 
log g= AH-B reffH-C rJf and providing in input the coefficients 
A, B, C. This option is useful when the investigated stars be- 
long for instance to the same stellar cluster and log g and Teff 
can be parametrized by a simple relation (i.e. as that described 
by a theoretical isochrone for a given evolutionary stage). 

The user can choose the way to treat log g (fixed or opti- 
mized following one of the methods described above); if the 
optimization of log g from the iron lines is requested but no 
Fe II lines are available, GALA will try to use the second op- 
tion (lines of other elements in different stages of ionization), 
or eventually will fix log g to the input value. 

5.2. Model atmospheres 

The algorithm used in GALA is basically independent 
from the code adopted to derive the abundances and from the 
model atmospheres. GALA is designed to manage both the 
two most used and publicly available models atmospheres, 
namely ATLAS9 and MARCS: 

• ATLAS9: The suite of Kumcz codes represents the 
only suite of open-source and free programs to face 
the different aspects of the chemical analysis (model 
atmospheres, abundance calculations, spectral synthe- 
sis), allowing any user to compute new models and up- 
grade parts of the codes. GALA includes a dynamic 
call to the ATLAS9 code 0. Any time GALA needs 
to investigate a given set of atmospheric parameters, 
ATLAS9 is called, a new model atmosphere is com- 
puted and finally it is stored in a directory. The latter is 
checked by GALA whenever a model atmosphere is re- 
quested, and ATLAS9 called only if the model is lack- 
ing. The convergence of the new model atmosphere 
is checked for e ach atmospher ic layer Following the 
prescriptions by ICastefiH (119881) . we require errors less 
than 1% and 10% for the flux and the flux derivative, 
respectively. Additional information about the calcula- 
tion for each model atmosphere is saved. In the current 
version, G ALA is able to manage th e grid of ATLAS9 
models by ICastelh & Kuruc j (120041) and the new grid 

^ T he ATLAS9 source code is available at the website 
|http://wwwuser.oat. ts.astro .it/castelli/sources/atlas9cod es. html I 



of models calculated by iMeszaros et al.l (120 12h for the 
APOGEE survey. 

• MARCS: At variance with ATLAS9, for the MARCS 
models the code to compute new model atmospheres 
is not released to the community. However, the Upp- 
sala group provides a large grid of the MARCS models 
on their websitfl When GALA works with these grids 
(including both plane-parallel and spherical symmetry), 
new models are computed by interpolating into the Up- 
psala grid by us ing the code developed by T. Masseron 
(lMasseron! '2006) 0. This code has been modified in or- 
der to put the interpolated MARCS models in ATLAS9 
format to use with WIDTH9. 

Note that the automatization of the chemical analysis based 
on EWs needs a wide grid of model atmospheres (both to in- 
terpolate and compute new models) linked to the code in order 
to freely explore the parameter space. Thus, GALA is linked 
to the ATLAS9 grid both with solar-scaled and a-enhanced 
chemical composition and to the MARCS grid with standard 
composition. Also, the use of other models or model grids 
can be easily implemented in the code. Sometimes, peculiar 
analysis or tests need to use specific models, for instance for 
the Sun (see the set of solar model atmospheres available in 
the website of F. Castelli) or with arbitrary chemical compo- 
sitions, as those computed with the ATLAS 12 code ( Castellil 
2005b). When a specific, single model is called, all the opti- 
mization options are automatically swichted off. 

5.3. Exit options 

GALA is equipped with a number of exit flags in order to 
avoid infinite loops or unforeseen cases stopping the analysis 
of the entire input list of stars. We summarize here the main 
exit options: 

( 1 ) the user can set among the input parameters the maximum 
number of allowed iterations for each star When the code 
reaches this value it stops the analysis, moving to the next 
star Generally this parameter depends on the adopted grid 
steps and if the input atmospheric parameters are close or not 
to the real parameters. When the Guess Working-Block is 
used, typically the analysis block converges in 3-5 iterations; 

(2) if the dispersion around the mean of the abundances of 
the lines used for the optimization (after the line rejection) 
exceeds a threshold value GALA skips the star In fact, very 
large dispersions can possibly suggests some problems in the 
EW measurements; 

(3) if the number of the lines used for the optimization (after 
the line rejection) is smaller than a threshold, the optimization 
is not performed and the atmospheric parameters are fixed to 
the input guess values; 

(4) the procedure is stopped if the requested atmospheric pa- 
rameter is outside the adopted grid of model atmospheres; 

(5) GALA skips the analysis of the star if the call to the model 
atmosphere fails (problems in the ATLAS9 models computa- 
tion or in the MARCS models interpolation) and the required 
model is not created. Otherwise, if the ATLAS9 model is 
calculated but some atmospheric layers do not converge (ac- 
cording to the criteria discussed above) GALA continues the 
analysis but it advises the user of the number of unconverged 
layers. 

^ http://marcs.astro.uu.se/ 
The original code is available at |http ://marcs .astro.uu. se/ software .php | 
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6. UNCERTAINTIES 

The code is equipped with different recipes to compute the 
uncertainties on each derived abundance. Several sources of 
error can affect the determination of chemical abundances, 
mainly the uncertainties due to the EW measurements and to 
the adopted \oggf (that are random errors from line to line) 
and those arising from the choice of the atmospheric param- 
eters (that are random errors from star to star but systematic 
from line to line in a given star). These uncertainties are quan- 
tified by GALA while other sources of errors (as the choice of 
the abundance calculation code or the adoption of the grid of 
model atmospheres) are neglected because considered exter- 
nal errors. 

6.1. Abundances statistical errors 

The statistical uncertainty on the abundance of each ele- 
ment is computed by considering only the surviving lines af- 
ter the rejection process (see Sect. 2)- When the uncertainty 
on the EW is provided for each individual line (thus allowing 
to compute the abundance error for each transition), the mean 
abundance is computed by weighing the abundance of each 
line on its error, otherwise simple average and dispersion are 
computed. For those elements for which only one line is avail- 
able, the error in abundance is obtained by varying the EW of 
laEw (if the EW uncertainties are provided). Otherwise, the 
adopted value is zero. As customary, the final statistical error 
on the abundance ratios is defined as <J j ^JNimes- 

6.2. Uncertainties on the atmospheric parameters 

GALA estimates the internal error for each stellar parame- 
ter that has been derived from the spectroscopic analysis. The 
uncertainties of Teff, Vt and log g are estimated propagating the 
errors of the corresponding optimization parameter: 



CTXi 



^C(Xi) 

^ SX; > 



where are T^ff, log g and v,, C(Xi) indicates the optimiza- 
tion parameters defined in Section ITTI and (Jc{Xi) the cor- 
responding uncertainties. 

The terms ^f^'^ (which parametrize how the slopes and the 
iron difference vary with the appropriate parameters) are c al- 
culated numerically, by varying Xi locally around the final 
best value, assuming the step used in the optimization process 
and recomputing the corresponding C(Xi). 

The terms C(Xi) are co mputed by app lying a JackJmife 
bootstrapping technique (see lLuptonll993h . The quoted quan- 
tities are recomputed by leaving out from the sample each 
time one different spectral line (thus, given a sample of N 
lines, each C(X) is computed N times by considering a sub- 
sample of N-1 lines). The uncertainty on the parameter X 
is (Tjack = vN-loSub, where asub is the standard deviation of 
the C(X) distribution derived from N sub-samples. The ajack 
takes into account the uncertainty arising from the sample size 
and the line distribution, and this resampling method is espe- 
cially useful to estimate the bias arising from the lines statis- 
tics. Note that the computation of the slopes is performed 
taking into account the effect of the EW and abundance un- 
certainties of each individual line. Thus, the uncertainty in the 
atmospheric parameter X; becomes 



It is worth to notice that these uncertainties represent the 
internal error in the derived parameters and are strongly de - 
pendent on the number of used lines and on the distribution 
of the lines (weak and strong transitions for the estimate of 
Vt and low and high x lines for Teff). Other factors that can 
affect the determination of the parameters (for instance, the 
threshold adopted in the EWs and in (Tew) are not included in 
the error budgets and they can be considered as external er- 
rors. Finally, the error due to the adopted grid size could be 
considered as a systematic uncertainty (being the same for all 
the analysed stars) and eventually added in quadrature to the 
internal error estimated by GALA. 

6.3. Abundances uncertainties due to the atmospheric 
parameters 

The evaluation of the uncertainties arising from the atmo- 
spheric parameters is a more complex task. Generally these 
errors are referred to as "systematic" uncertainties but this 
nomenclature is rather imprecise. In fact, the variation of a 
given parameter changes the abundance derived from the lines 
of the same element in a similar way (for instance, an increase 
of Tgs increases the abundance of all the iron lines). However, 
this error will be different from star to star, due to the different 
number of lines, strength and x distributions, EWs quality and 
so on. Thus, the uncertainties from the atmospheric parame- 
ters should be considered as random errors, when we compare 
different stars (but they are systematic uncertainties from line 
to line). 

Several recipes are proposed in the literature. The most 
common method is to re-compute the abundances changing 
each time one parameter only, and keeping fixed the other 
ones to their best estimates. Then, the corresponding vari- 
ations in the abundances are added in quadrature. This ap- 
proach is the most conservative, because it neglects the co- 
variance term s ari sing from the interplay among the parame- 
ters (see Sect. \2.2\ . providing only an upper limit for the total 
error budget. 

GAL A follows the approach described by ICayrel et al.l 
(120041) to naturally take into account the covariance terms. 
When the optimization process is ended, the analysis is re- 
peated by altering the final Teff by h-ctt;,, and -aj^^ (these un- 
certainties are calculated as described in Section 6.2), and re- 
optimizing the other parameters. The net variation of each 
chemical abundance with respect to the original value is as- 
sumed as final uncertainty due to the atmospheric parame- 
ters and including naturally the covariance terms. Addition- 
ally, under request, also the abundance variations following 
the classical approach to vary one only parameter each time 
(keeping the other parameters fixed) are calculated, leaving 
the user free to use this information as preferred. 

6.4. Quality parameter for the final solution 

GALA provides also a check parameter, useful to judge the 
quality of the global solution and to identify rapidly stars with 
unsatisfactory solutions. For each model used during the op- 
timization process a merit function Fmcnf is defined as: 
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taking into account the values of the optimization parame- 
ters and the corresponding uncertainties. In an ideal case, 
F„,erit is zero if the three optimization parameters are exactly 
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zero. Generally, all the solutions with Fmerit — 1 are valid and 
equally acceptable, while values of Fmerit > > 1 are suspect and 
point out that at least one of the parameters is not well con- 
strained within the quoted uncertainty. Note that Fmerit pro- 
vides only an indication if the solution is acceptable or not, 
but it does not specify which parameter is not well defined. 

Summaryzing, when the full optimization process is com- 
pleted, GALA will provide for each analysed element the 
(weighted) mean abundance, the dispersion and the number 
of used lines (which provide the statistical uncertainty), the 
net variation in abundance due to the new optimization with 
Teff + cr%„ that with Teff - (Tj;,, (which provide the uncer- 
tainty owing to the choice of stellar parameters). Also, for 
each atmospheric parameters the quoted internal uncertainties 
are computed. Finally, the quality parameter /Vnenv is provided 
to evaluate the goodness of the solution as whole. 

7. DEPENDENCE ON SNR 

We performed a number of experiments to test the stabil- 
ity and reliability of the derived atmospheric parameters with 
GALA at different noise conditions. We performed two kind 
of experiments, described in the following: the first based on 
a grid of synthetic spectra of abundances and atmospheric pa- 
rameters known a priori, in order to estimate the reliability 
of the code as a function of the parameters and the signal- 
to-noise; the second group of tests is based on real spectra 
already analysed in literature. In the foll owing, the EWs were 
measured by means of DAOSPEC (Stet son & Pancinoll2008h 
adopting a Gaussian profile for the line fitting. 

7.1. Synthetic spectra at different noise conditions 

We analysed with GALA a grid of synthetic spectra, com- 
puted to mimic the UVES@VLT high resolution spectra with 
the 580 Red Arm setup. The grid of synthetic spectra includes 
SNR of 20, 30, 50, 100 per pixel for two different sets of at- 
mospheric parameters: reff= 4500 K, logg=1.5, Vt= 2 km/s, 
[M/H]= -1.0 dex to simulate a giant star, and Tef{= 6000 K, 
logg=4.5, Vt= 1 km/s, [M/H]= -1.0 dex to simulate a dwarf 
star The spectra were computed with the following proce- 
dure: 

(1) for a given model atmosphere, two synthetic spectra were 
calculated with the SYNTHE code over the wavelength range 
covered by the two CCDs of the 580 UVES Red Arm grating 
and then convoluted with a Gaussian profile in order to mimic 
the formal UVES instrumental broadening; 

(2) the spectra were rebinned to a constant pixel-size 
(6X= 0.0147 and 0.0174 pixel/A for the lower and upper chips 
respectively); 

(3) the synthetic spectra (normalized to unity) were multiplied 
with the efficiency curve computed by the FLAMES-UVES 
ESO Time Calculator in order to model the shape of the tem- 
plates as realistically as possible; 

(4) Poissonian noise was injected in the spectra to simulate 
different noise conditions. Basically, the SNR varies along 
the spectrum, as a function of the efficiency (and thus of the 
wavelength). The noise was added in any spectrum accord- 
ing to the curve of SNR as a function of A provided by the 
FLAMES-UVES ESO Time Calculator For each SNR a sam- 
ple of 200 synthetic spectra were generated. 

Fig. |8] summarizes the average values obtained for each 
MonteCarlo sample for each atmospheric parameter as a func- 
tion of SNR; the errorbars indicate the dispersion around the 
mean. Results of the simulations of the giant star model atmo- 
sphere are shown in the upper panels of each window, while 



the lower panels summaryze the results for the dwarf star sim- 
ulations. Basically, the original parameters of the synthetic 
spectra (marked in Fig. |8]as dashed horizontal lines) are re- 
covered with small dispersions and without any significant 
bias. The major departure from the original values is ob- 
served in the microturbulent velocities (both dwarf and giant) 
at SNR= 20, because of the loss of weak Unes. 
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Fig. 8. — Average values for the recovered atmospheric parameters as a 
function of SNR for the MonteCarlo samples described in Sect. |7.1| upper 
panels of each window show the results for the giant star model, while lower 
panels show the results for the dwarf" stars. Errorbars are the dispersion by 
the mean. Dashed hnes are the original value for each parameter 



7.2. Real echelle spectra at different noise conditions 

The previous test provides an indication about the stability 
of the method against the SNR, starting from spectra their pa- 
rameters are known a priori. However, the employed noise 
model is a simplification because it does not take into account 
some effects that can also heavily affect the measurement of 
the EWs, such as the correlation of the noise among adjacent 
pixels, flat-fielding residuals, failures in the echelle orders 
merging, presence of spectral impurities. Also, the atomic 
data of the analysed lines are the same used in the compu- 
tation of the synthetic spectra, thus excluding from the final 
line-to-line dispersion the random error due to the uncertainty 
on the atomic data. 

In order to provide an additional test about the perfor- 
mance of GALA in conditions of different noise, we per- 
formed a simple experiment on the spectra acquired with 
UVES@FLAMES of the giants stai- NGC 1786-1501 in the 
Large Magellanic Clo ud globular cluster NGC 1786 (see 
iMucciarelli etal]|2009l 12010). This is a sample of 8 spectra 
with the same exposure time (^45 min) obtained under the 
same seeing conditions was secured, with a typical SNR per 
pixel of 20 for each exposure. We used this dataset in order to 
obtain 8 spectra with different SNR ranging from ^20 to ~'60 
depending on the number of exposures averaged: the spec- 
trum with the lowest SNR is just one exposure, the spectrum 
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with the largest SNR the average of all the 8 acquired expo- 
sures. The derived parameters for each spectrum are showed 
in Fig. |9] as a function of SNR. The errorbars are derived by 
applying the Jackknife bootstrap technique for T^ff, log g and 
Vt, while for [Fe/H] we used the dispersion by the mean as 
estimate of the error 

We note that the parameters are well constrained with small 
uncertainties for spectra with low SNR: this is not a numer- 
ical artifact of the code but it is due to the large number of 
transitions available in the UVES spectra, coupled with an ac- 
curate rejection of the outliers and the use of the uncertainties 
for each individual line in the slopes computations. Also, we 
note that major departures from the final parameters are again 
in the determination of Vt at SNR= 20: the derived low value 
of Vt is due to the fact that at low SNR several weak lines 
(useful to constrain the microturbulent velocity) are not well 
measured (then discarded by GALA) or not identified in the 
noise envelope by DAOSPEC. Note that the derived trend of 
Vt as a functi on of SNR shows t he same behaviour found and 
discussed bv lMucciarellil (1201 lb . 
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8.1. Selection of the lines 

For each star we defined a suitable linelist of Fe 1 and Fe II 
transitions, starting from the most updated version of the Ku- 
rucz/Castelli lines dataset^]- We apply an iterative procedure 
to define the linelist. Assuming that the parameters of the 
targets are not known a priori, we performed a first analysis 
by using a preliminary linelist including only laboratory 
transitions with x <6 eV and log (gf)>-5 dex. Such a linelist 
is not checked against the spectral blendings arising from the 
adopted spectral resolution and the atmospheric parameters 
and it is used only to perform a preliminary analysis. With the 
derived new parameters we define a new linelist for each star. 
The lines are selected by the inspection of synthetic spectra 
computed with the new parameters and convoluted with a 
Gaussian profile in order to reproduce the observed spectral 
resolution. At this step, only iron transitions predicted to 
be unblended are taken into account and used for the new 
analysis. 



8.2. EWs measurements 

EWs are measured by using the code DAOSPEC which 
adopts a saturated Gaussian function to fit the line profile and 
an unique value for the full width half maximum (FWHM) 
for all the lines. We start from the FWHM derived from the 
nominal spectral resolution of the spectra, leaving DAOSPEC 
free to re-adjust the value of FWHM according to the global 
residual of the fitting procedure. The measurement of EWs is 
repeated by using the optmized FWHM value as a new input 
value, until convergence is reached at a level of 0. 1 pixel. The 
formal error of the fit provided by DAOSPEC is used as 1 cr 
uncertainty on the EW measurement. 
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Fig. 9. — Behaviour of the derived atmospheric pai'ameters for the giant star 
NGC 1786-1501 as a function of SNR obtained by using different UVES co- 
added spectra. Errorbars are derived from the Jackknife bootstrap technique 
for Teff, log g and I't, and as dispersion by the mean for [Fe I/H]. 



8. AN EFFICIENT APPROACH: ARCTURUS, SUN, HD 84937 AND 
/xLEONIS 

In this section we describe a convenient and robust method 
for performing abundance analysis with GALA, applied to the 
case of 4 stars (namely, the Sun, Arcturus, HD 84937 and 
/iLeonis) of different metallicity and evolutionary stage and 
whose parameters are well established among the closest F-G- 
K stars. We retrieved high-resolution (^45000) spectra from 
the ESC0 (for Sun, Arctums and HD 84937) and ELODIEH 
(for //Leonis) archives. 

' http://archive.eso.org/eso/eso_archive_main.html 
http://atlas.obs-hp.fr/eIodie/ 



8.3. Analysis with GALA 

The programme stars are analysed by employing all the 
three Working-Blocks, requiring a spectroscopic optimization 
of Teff, log g, Vt and [M/H] and starting in all cases from the 
same set of guess parameters (namely reff= 5000, logg= 2.5, 
[M/H]= -1.0 dex and Vt= 1.5 km/s). The optimization is per- 
formed by exploring the parameters space in small steps of 
5Tef{= 50 K, Slog g= 0.1 and Svi= 0.1 km/s; the metallic- 
ity is investigated by adopting the step of the ATLAS9 grids 
((5[M/H]= 0.5 dex). 

In the first run, we analysed the programme stars assuming 
the same configuration for the input parameters of GALA, in 
particular we included only lines with aEw <10% and with 
EWR>-5.8 (coiTesponding to --10 mi at 6000 A). After a 
first run of GALA, we refined the maximum allowed EWR, 
that depends mainly by the onset of the saturation along the 
curve of growth (and thus it is different for stars with different 
atmospheric parameters). We adopted as maximum allowed 
value, EWR= ^.65 for Arcturus and /iLeonis, and ^.95 for 
the Sun and HD 84937. These values are chosen on the basis 
of the visual inspection of the curve of growth, in order to ex- 
clude too strong lines, for which the Gaussian approximation 
can fail. After a first run of GALA, the linelist is refined by us- 
ing the new parameters obtained by GALA as described above 
and the procedure repeated. Table 1 summaryzes the de- 
rived atmospheric parameters (with the corresponding Jack- 
knife uncertainties) and the [Fe/H] I and [Fe/H] II abundance 

' ' http://wwwuseroat.ts.astro.it/castelli/linelists.html 
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ratios, together with the two errorbars due to the atmospheric 
parameters and the internal error computed as cr / ^/iNnnes)- 

We compare our results with those available in literature 
(and listed in Table 1 as reference). We derive for the Sun 
7;ff= 5800±64 K, log g= 4.50±0.18, \\= 1.20±0.13 km/s 
and [Fe/H]=-0.01±0.03 dex (where the errorbar is the sum 
in quadrature of the individual uncertainties listed in Table 
1). Our results for Teff and log g agree very well with those 
listed in the compilation of the NASA websiteQ Concern- 
ing the microtubulent velocities, values available in litera- 
ture range fro m 0.8 km/s (Biemontet al. 1981) to 1.35 km/s 
jSteffen.Ludwi'g & C affau 2009). A value of 1 km/s is typi- 
cally adopte d as representative f or the Sun in several chemical 
analysis (see lCaffau et al.ll20Tlh . 

Our analysis of Arcturus provides Tsif= 4300±60 K, 
log g=1.60±0.06 and Vt= 1.50±0.06 km/s, with an iron abun- 
dance [Fe/H]= -0.51±0.05. These results w ell match the re- 
cent analysis of Arcturus by Ramirez & Allende Prietol (l201 lb 
that provide an accurate determination of the atmospheric pa- 
rameters and the chemical composition; in particular T^s and 
log g are derived in an independent way with respect to our 
approach, finding reff= 4286±30 K (by fitting the observed 
spectral energy distribution), log g=1.66±0.05 (through the 
trigonometric parallax), while Vt turns out to be 1 .74 km/s (by 
using the same approach used in GALA). The final iron abun- 
dance is [Fe/H]= -0.52±0.04 dex. In both cases, the agree- 
ment with the literature values is good. Finally, Fig.fTOlshows 
as example the graphical output of GALA for Arcturus. 



s? 7.5 

OJ 

-d 

^ 7 

0) 

V 6,5 
6 

6? 10 



I I I I I I I I I I I I I I 



-5.5 -5 
log(EW/A) 





5500 6000 6500 
A(A) 



< -5 



2 4 
x(eV) 



\ 

h 

-7 -6 -5 

iog{gf)-ex 



Fig. 10. — Example of the graphical output of GALA for Arcturus (see 
Section[8): black circles are the Fe I lines used in the analysis and the empty 
circles show rejected points. 

For HD 84937 we derive %ff= 6150±56 K, 
log g= 3.20±0.13, vt= 0.70±0.24 km/s and [Fe/H]= - 
2.28±0.06 dex, while for /iLeonis we obtain reff= 4500±81 
K, log g=2.40±0.26, vt= 1.40±0.07 km/s and 
[Fe/H]= H-0.37±0.06 dex. For these two stars several 

http://nssdc.gsfc.nasa.gov/planetary/factsheet/sunfact.html 



determinations are available in literature and we decide to use 
as reference the average of the values hsted in the classica l 
compilation by ^ Cavrel de Strobel. Soubiran & Ralitel (12001b . 
providing 7;ff= 6251±94 K, log g= 3.97±0.18 and [Fe/H]= - 
2.14±0.17 dex for HD 84937 and T^s= 4504±121 K, 
log g= 2.33±0.27 and [Fe/H]= H-0.28±0.13 dex for ^Leonis. 
No te that the value of the microturbulent veloc ities is omitted 
by ICavrel de Strobel. Soubi ran & RaUtel (120011) because the 
listed authors in their compilation use different definitions 
for this parameter (see 'Muccia relliil201 1[ for a review of the 
different approaches) or assume a representative value. For 
HD 84937 the values of Vt range from 0.8 up to 1.7 km/s, 
while our value is slightly lower Also for /xLeonis the range 
of values for Vt is wide (from 1 .2 up to 2.2 km/s) but our value 
agrees with this value. Basically, the agreement with the 
values listed by Cayrel de Strobel. S oubiran & Ralite (200 ij) 
is good, but for the gravity of HD 84937, for which we derive 
a lower value; this difference can be partially explained in 
light of the different Vt. 

8.4. Stability against the initial parameters 

A relevant feature of an automatic procedure to infer param- 
eters and abundances is its stability against the input atmo- 
spheric parameters. In order to assess the effect of different 
first guess parameters, we analyse the spectrum of Arcturus 
by investigating a regular grid of input parameters with Teff 
ranging from 3800 to 4800 K (in steps of 100 K) and log g 
from 1.0 to 2.2 (in steps of 0.1). Fig.fTTIshows the grid of the 
input parameters in the Teff-log g plane (empty points) with 
the position of the derived parameters (black points); the up- 
per panel summarizes the results when GALA is used without 
the Refinement Working-Block, while the lower panel shows 
the results obtained by employing also the refinement option. 
The recovered parameters cover a small range: in the first run 
the dispersion of the mean is of 46 K for T^ff and 0.09 for 
log g, while these values drop to 25 K and 0.05 respectively, 
when the Refinement Working-Block is enabled. 

9. A TEST ON GLOBULAR CLUSTERS 

Globular clusters are ideal templates to check the capability 
of our procedure deriving reliable atmospherical parameters, 
because of the homogeneity (in terms of metallicity, age and 
distance) of their stellar content. Thus, the derived parameters 
for stars in a given globular cluster can be easy compared with 
theoretical isochrones in the Teff-log g plane. 

We apply the same procedure described in Sect.|8]to analyse 
a set of high-resolution spectra for stars in the globular cluster 
NGC 6752, ranging from the Turn-Off up to the bright por- 
tion of the Red Giant Branch. The spectra have been retrieved 
by the ESO archiv4]3 and reduced with the standard ESO 
pipelin43- They are from different observing programmes 
and with different SNR, including giant stars crossing the Red 
Giant Branch Bump region observed with UVES@VLT (slit 
mode) within the ESO Large Program 65.L-0165 with very 
high (>200) SNR, the stars in the bright portion of the Red 
Giant Branch observed with UVES -FLAMES @VLT (fiber 
mode) within the Galactic globular clusters survey presented 
by Camtta et al. (2009) and the dwarf/subgiant stars observed 
with UVES@VLT (slit mode) within the ESO Large Program 
165.L-0263. 

" http://archive.eso.org/cms/eso-data.html 
http://www.eso.Org/sci//software/pipelines/ 
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TABLE 1 

Iron abundance (from neutral and singly ionized lines) and atmospheric 

PARAMETERS DERIVED WITH GALA FOR ARCTURUS, /^Leonis, THE SUN AND HD 84937. 
For the ABUNDANCES, THE FIRST TWO ERRORBARS ARE THE UNCERTAI NTIES ARISING 
FROM THE ATMOSPHERIC PARAMETERS FOLLOWING THE PRESCRIPTIONS BY'CAYREL ET ALJ 
(3004;), WHILE THE LAST ONE IS THE INTERNAL ERROR CALCULATED AS cr/ v'(A';„,„). THE 
UNCERTAINTIES IN THE DERIVED ATMOSPHERIC PARAMETERS ARE COMPUTED WITH A 
JACKKNIFE bootstrap TECHNIQUE. IN THE LOWER PART OF THE TABLE, THE VALUES 
AVAILABLE IN LITERATURE ARE LISTED FOR COMPARISON. 



Star 


[Fe I/H] 
(dex) 


[Fe II/H] 
(dex) 


(K) 


logg 


Vl 

(km/s) 




Sun 
Arcturas 
HD 84937 
^Leonis 


-0.0l!i;"^± 0.009 
-0.5 l^H'g'i 0.007 
-2.28^-'|6 ± 0.007 
+0.37S± 0-011 


+0.01^ "1 ±0.003 
-0.52™-|;3 ±0.015 
-2.27:t<'-''5 ±0.020 
+0.38™ ±0.033 


5800±64 
4300±60 
6150±56 
4500±81 


4.50±0.18 
1.60±0.06 
3.20±0.13 
2.40±0.26 


1.20±0.13 
1.50±0.06 
0.70±0.24 
1.40±0.07 


Sun 
Arcturas 
HD 84937 
/^Leonis 


0.00 
-0.52±0.02 
-2.14±0.17 
+0.28±0.13 


0.00 
-0.40±0.03 


5778 
4286±30 
6251 ±94 
4504±121 


4.438 
1.66±0.05 
3.97±0.18 
2.33±0.27 
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Fig. 1 1 . — Position of the final parameters for Arctums (black points) in 
the Teff-logg plane, in comparison with the input parameters (empty circles), 
obtained by using GALA without (upper panel) and with (lower panel) the 
Refinement Working-Block (upper and lower panel, respectively). 

The main panel of Fig. [12] shows the position of the fi- 
nal parameters derived with GALA in the Teff-log g plane. 
Also, two theoretical isochrones with an age of 12 Gyr and a 
metallicity of Z= 0.0006 (assuming an a-enhanced chemical 
mixture) are shown as reference (gre y curve is from BaSTl 
database by iPietrinfe rni et al .l (120061) and black curve from 
Padua database by|Girardi et al.l l2000l) . In the lower panel, 
we show the behaviour of the [Fe/H] ratio as a function of 
Teff. No significant trend is found, while the star-to-star scat- 
ter increases with increasing T^ff because of the lower SNR. 

The parameters derived with GALA well reproduce the be- 
haviour predicted by the theoretical models for a old simple 
stellar population with the same metallicity of the cluster, con- 
firming the physical reliability of the final solution. Also, we 
note that the errorbars, both in Teff and log g, change according 
to the quaUty of the spectra, ranging from ~30 K and ~0.15 



; NGC 6752 

1 - BaSTI Padua 

_ Age = 12 Gyr, Z = 0.0006 




5000 
T.„ (K) 



Fig. 12. — Main panel: position in the Teff-log g plane of the stars in 
the globular cluster NGC 6752 analysed with GALA . We plotted as refer- 
ences two isochrones computed with an age of 12 Gyr and a metallicity of 
Z= 0.0006 (assuming an Q-enhancement ctiemical mixture), from the BaSTI 
i Pietrinfemi et al. 2006, dotted curve) and Padua (Girardi et al. 2000, solid 
curve) database. Lower panel: behaviour of the [Fe/H] ratio as a function of 

Tcff. 

for the giants with the highest SNR up to ~200 K and --0.5 
for the dwarf stars with the lower SNR. 

The entire sample of 52 stars provides an average iron 
abundance of [Fe/H]l=-1.63 dex (a= 0.06 dex). This value is 
consistent with the previous estimates available in literature 
that point out a n iron content ranging from [Fe/H] =-1.62 
jGrundahl et"ani2002h up to [Fe/H]=-1.42 dex dGratton et al] 
200l!)- In this comparison we cannot take into account the 

different adopted solar values. 

The stars in common with Carretta et al.l (120091) show a 
reasonable agreement in the atmospheric parameters, with 
T^fj^ -T^f f""'=+A9 K (cr= 42 K), logg^'^^-logg^'"'""''=- 

0.26 (a= 0.07) and [Fe/H]^^^-[Fe/Hf'"''''"'"=-Om dex 
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(cr= 0.05 dex). 

Also for the stars in common with lYong et al] (l2005h the 
agreement is excellent (also because these spectra have very 
high SNR and a large wavelength coverage, thus allowing 
the measurements of a large number of Fe 1 and Fe 11 lines): 

T^'jj^-Tj/^>^=-\ K (cr= 32 K), logg'^^^ -I ogg'^""«=-0. 39 

(cr= 0.11) and [Fe/H]'^^^-[Fe/Hf°"s=+0.00 dex ia= 0.06 
dex). 

The comparison with the a nalysis of the Turn-Off and 
Sub-Giant Branch stars by iGratton et all (12001b is not 
trivial, because they derived the atmospheric parameters 
from median spectra for the two groups of stars. The 
temperatures have been derived by fitting the wings of 
the Ha, while the gravities have been obtained from the 
position of the stars in the Color- Magnitude Diagram . 
When we compared the parameters by IGratton et al.l (12001b 
with the average values obtained by our analysis of indi- 
vidual stars, the agreement is not perfect but consistent 
within the uncertainties: we find < T^'jj-^ >-T^^/f"""=+l64 

K, < logg'^^^ >-logg'^'''"'°"=+0.06 and [Fe/Hf"^^- 
[Fe/Hf ''""'"=-0.19 dex (cr= 0.15 dex) for the Turn-Off stars, 
and < T^'jj^ >-T^/}"""'=-\5 K, < logg^'^^ >-logg^''""""=- 
0.36 and {FelH]^'^^-{Fe/H]^''"'""=-Q.21, dex {a= 0.11 dex) 
for the Sub-Giant stars. 

10. SUMMARY 

In this paper we have presented a new, automatic tool to per- 
form accurate analysis of stellar spectra. GALA is designed 
to perform automatically the search for the best atmospheric 
parameters Teff, log g, V/ and the overall metallicity [M/H] 
for moderate and high resolution stellar absorption spectra, 
by using the EWs of metallic lines. Also, GALA provides the 
abundance of each individual line for which the user provides 
the EW, as well as the average abundance for each atomic 
species. 

The source code of GALA is freely available at the web- 
site http://www.cosmic-lab.eu/Cosmic-Lab/Products.html to- 
gether with the user manual (including the information about 
installation, configuration of the input files and how to obtain 
and properly use the model atmospheres) and an example of 
the input files as reference. The main advantages of the code 
are the capability: 

1 . to optimize all the parameters or only part of them. The 
code is versatile in order to perform different kind of 
analysis (full or partial spectroscopic analysis, experi- 
ments about the guess parameters...) and adopting dif- 
ferent recipes to derive log g; 

2. to perform a careful rejection of the outliers according 
to the line strength, the EW quality and the line distri- 
bution in the A(Fe)-x and A(Fe)-EWR planes; 

3. to estimate for each individual star the internal errors 
for (a) the optimized parameters by adopting the Jack- 
knife bootstrapping technique, and (b) the derived un- 
certainties due to the choice of atm ospheric parameters , 
following both the prescriptions bv lCavrel et al.l (12004b 
and the classical method of altering one parameter at a 
time. 

We have performed an extensive set of test with both syn- 
thetic and observed spectra, in order to assess the perfor- 
mances of the code. 



1 . Experiments with synthetic spectra (whose atmospheric 
parameters are known a priori) with the injection of 
Poission noise to simulate different noise conditions 
show a high stability of the code to recover the atmo- 
spheric parameters, without significant bias. The major 
departure from the original values is found in the mi- 
croturbulent velocity of low (~^20) SNR spectra, due to 
the systematic loss of weak lines. 

2. A set of FLAMES-UVES spectra of the LMC giant 
star NGC 1786-1501 observed with different signal-to- 
noise (from ^20 up to ^60) has been analysed with 
GALA, confirming that our procedure well constrains 
the parameters also in case of low spectral quality. Also 
in this case, we found that the largest departures are 
for the microturbulent velocity in the spectra with low 
SNR, because the weak lines are not well measured or 
not detectable in the noise envelope, leading to an un- 
derestimate of this parameter. 

3. We analysed 4 stars (namely, the Sun, Arcturus, 
HD 84937 and /xLeonis) of different metallicity and 
evolutionary stage and whose parameters are well es- 
tablished among the closest F-G-K stars. We described 
an efficient method (including the line selection, the 
measurement of the EWs and the chemical analysis) to 
best exploit the capabilities of GALA. Our results for 
these stars (both for atmospheric parameters and [Fe/H] 
ratio) well agree with those available in literature. 

4. Finally, we analysed a sample of 52 stars of the Galac- 
tic globular cluster NGC 6752, in different evolutionary 
stages. The derived Tgjf and log g well follow the pre- 
dictions of theoretical isochrones with the appropriate 
age and chemical composition for this cluster Also, 
the derived iron content nicely agrees with the previous 
estimates of other works. 

The code permits to obtain chemical abundances and atmo- 
spheric parameters for large stellar samples in a very short 
time, thus making GALA an useful tool in the epoch of the 
multi-object spectrographs and large surveys. Because of its 
nature of open source code, GALA will be implemented in 
the next releases according to the feedback with the users. In 
particular, we plan to include in the code new grids of ODFs 
and models that will be publicly released in the future. Also, 
GALA will be constantly updated in order to include varia- 
tions and changes in ATLAS9 and MARCS models, as well 
as in the ATLAS9 code. 

The authors warmly thank Fiorella Castelli, Thomas 
Masseron, Piercarlo Bonifacio and Andrea Negri for useful 
comments, discussions and suggestions. We thanks the 
anonymous referee for a careful reading of the paper and 
helpful comments. This research is part of the project 
COSMIC-LAB funded by the European Research Council 
(under contract ERC-2010-AdG-267675). 

Und wenn dich das Irdische vergafi, zu der stillen Erde sag: 
Ich rinne. Zu dem raschen Wasser sprich: Ich bin. (R. M. 
Rilke) 
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